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Using density functional theory calculations including a Hubbard U term we explore the effect of strain and 
confinement on the electronic ground state of superlattices containing the band insulator LaAlOs and the cor- 
related metal LaNiO:?. Besides a suppression of holes at the apical oxygen, a central feature is the asymmetric 
response to strain in single unit cell superlattices: For tensile strain a band gap opens due to charge dispro- 
portionation at the Ni sites with two distinct magnetic moments of 1.45/iB and 0.71/ib- Under compressive 
stain, charge disproportionation is nearly quenched and the band gap collapses due to overlap of d 3z 2_ r 2 bands 
through a semimetallic state. This asymmetry in the electronic behavior is associated with the difference in oc- 
tahedral distortions and rotations under tensile and compressive strain. The ligand hole density and the metallic 
state are quickly restored with increasing thickness of the (LaA103)„/(LaNi03)„ superlattice from n = 1 to 
n = 3. 



PACS numbers: 73.20.-r,68.65.Cd,71.30.+h,74.78.Fk 

The ability to grow transition metal oxide heterostructures 
with atomic control of layer thickness opens an exciting realm 
to explore confinement effects and novel effectively two- 
dimensional electronic phases that are not accessible in the 
bulk compounds. To this end, Chaloupka and Khaliullin 
recently proposed that a cuprate-like behavior may be stabi- 
lized in a La2NiA/06 superlattice (SL) where MO2 and MO2 
planes alternate along the c-axis of the ABO3 perovskite lat- 
tice and M denotes a trivalent cation (M=A1, Ga, Ti). LaNiC>3 
is a perovskite compound with rhombohedral symmetry (i?3c) 
that remains metallic at all temperatures due to degenerate e g 
orbitals at the Ni 3+ -sites (t% g e x ,S = 1/2). In contrast to bulk 
LNO, all other representatives of the rare earth nickelates un- 
dergo a metal-to-insulator transition (MIT) at low tempera- 
tures, driven by charge ordering (CO) rather than the Jahn- 
Teller effect. v Chaloupka and Khaliullin proposed that the 
orbital degeneracy in LNO/LMO superlattices can be lifted 
and a selective d 32 2_ r 2 or d x 2_ y 2 occupation at the Ni 3+ sites 
achieved by applying strain induced by the substrate. Subse- 
quent density functional theory (DFT) calculations in the local 
density approximation (LDA) found a two sheet Fermi surface 
(FS) with strongly hybridized d^ z 2_ r 2 and d x 2_ y 2 orbitals. ■ 
A single-sheet FS could only be obtained by including corre- 
lation effects within the dynamical mean field approximation 
(DMFT) and/or by using Sc as counterion to Ni instead of Al. 
Han et al. showed that the choice of the counterion (Al, Ga, 
B, In) rather than strain can be used to control chemically the 
orbital polarization. 

The DFT studies mentioned above have considered ex- 
clusively superlattices in a tetragonal setup. However, both 
LaNi0 3 (LNO) and LaA10 3 (LAO) are rhombohedral in the 
bulk (space group i?3c). Here we explore the effect of con- 
finement and strain on the octahedral distortions. DFT cal- 
culations with an on-site Coulomb repulsion term show that 
the lattice distortions have a significant influence on the band- 
width and thereby on the electronic ground state of the sys- 
tem. In particular, a charge disproportionated (CD) insulat- 
ing state is obtained in LAO1/LNO1 SLs under tensile strain. 
In a tetragonal cell such a state is reported only at higher U 



values. The tendency toward CD is vastly suppressed under 
compressive strain. 

A metal-to-insulator transition as a function of 
LNO thickness has been observed experimentally for 
LaNi0 3 /LaA10 3 1 " and LaNi0 3 /SrMn0 3 superlattices 
as well as in thin LNO films grown on an STO(OOl) sub- 
strate. • The Ruddlesden-Popper compound La4Ni 3 08 1415 
is another example for an insulating state due to confinement. 
Here, we explore the effect of quantum confinement by 
varying the period of the (LAO)„/(LNO)„ superlattice from 
n = 1 to n = 3. The results show a quick convergence toward 
bulk metallic behavior within the LNO part of the superlattice 
with increasing thickness. 



I. CALCULATIONAL DETAILS 

To gain insight into the electronic behavior of 
(LAO)„/(LNO)„ (n = 1,3) superlattices, we performed 
density functional theory calculations with the full potential 
linearized augmented plane wave method as implemented 
in the WIEN2k code. An on-site Coulomb repulsion term 
is considered beyond the generalized gradient approxi- 
mation of the exchange correlation potential within the 
LDA/GGA+[7 approach using U = 4 eV and J = 0.7 eV 
for Ni and U = 7 and J = eV for La. The muffin tin 
(MT) radii are 1.80 bohrs (Ni and Al), 2.30 bohrs (La) and 
1.60 bohrs (oxygen). Inside the muffin tins, wave functions 
are expanded in spherical harmonics up to = 10 and 
nonspherical contributions to the electron density and poten- 
tial are considered up to = 6. The energy cutoff for the 
plane wave representation in the interstitial is -E™^x = 25 Ry 
for the wave functions and E^ x = 144 Ry for the potential. 
76 fc -points were used for the integration in the irreducible 
part of the Brillouin zone (BZ). To model the effect of 
substrate- induced strain, the in-plane lattice parameter of 
the superlattice was set to that of LAO (a = 3.78A), STO 
(a = 3.905A) or DySc0 3 (DSO) (a = 3.94A), respectively. 
For the out-of -plane parameter we have used values obtained 
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from XRD measurements: 3.83 (STO) and 3.93 A (LAO). 
Relaxations of the c parameter for superlattices strained at 
a dso rendered a c value very similar to the one for the STO 
case. This implies a deviation from the volume conserving 
scenario for tensile strain, as also suggested from XRD 
measurements. While previous DFT studies ■ have used a 
tetragonal setup, we explored the influence of strain on the 
octahedral rotations by optimizing the internal parameters in 
a monoclinic unit cell containing 20 atoms for LNO bulk and 
the 1/1 superlattice. In the present study parallel alignment of 
the Ni spins is considered. Calculations for antiferromagnetic 
configurations were unstable and higher in energy. 
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FIG. 1: (Color online) Bulk LaNi0 3 . (a) Projected density of states 
of Ni 3d character (black dotted line) obtained within GGA+U 
(U — 4 eV and J = 0.7 eV). Positive (negative) values corre- 
spond to majority (minority) spin. Note that d 3z 2_ r 2 (green line) and 
d x 2 _ y 2 orbitals (orange dashed line) are degenerate in bulk LNO. (b) 
Electron density distribution of the unoccupied states integrated be- 
tween Ep and Ep+2.5 eV. La and O ions are denoted by purple and 
red spheres, respectively. 



Prior to studying the LAO/LNO superlattices we investi- 
gate here the electronic properties of LNO bulk. In order to 
facilitate comparison to the LNO/LAO superlattices the calcu- 
lations are performed in a monoclinic unit cell containing 20 
atoms. Figure 1(a) shows the projected density of states at the 
Ni sites. In bulk LNO d 3z 2_ r 2 and d x 2_ y 2 are degenerate and 
the broad e g band extends from -7 eV to 2 eV in the majority 
spin channel, while the minority e g states are not occupied ex- 
cept for some hybridization with the O 2p states between -4.5 
and -3.5 eV. The formal Ni valence in LNO is considered to be 
Ni 3+ (if s > e g)> but the occupation of the Ni e g states is much 
higher (both the d 3z 2_ r 2 and d x 2_ y 2 orbtails contain one elec- 
tron each; see Table I). As displayed in the hole density dis- 
tribution integrated between Ep and Ep + 2.5 eV [Fig 1(b)], 
the higher number of electrons at the Ni sites is counterbal- 
anced by holes on the ligands of p a character (~ 0.2e within 
the MT). This corresponds to a state closer to 3rf 8 L, where L 
denotes a hole on the oxygen site, as observed also for other 
rare-earth nickelates. For this reason bulk LaNi03 was ruled 
out as a possible parent compound for high-Tc superconduc- 



FIG. 2: (Color online) Band structure obtained within GGA+C7 of 
bulk LaNi0 3 along the high symmetry lines T(0, 0, 0) - Z(0, 0, it) - 
A(tt, 0, 7r) — R(tt, 7r, 7t) — Z(0, 0, 7r) using a \/2a x \/2a x 2a unit 
cell with a being the lattice parameter of the pseudocubic cell. In (a) 
and (b) the d 3z 2_ r 2 and d x 2_ y 2 character is emphasized by green 
triangles and light orange circles, respectively. 

The band structure of bulk LNO in Fig. 2 is plotted along 
the high-symmetry lines in the BZ of the 20-atom cell. We 
concentrate on the antibonding e g bands in the majority spin 
channel extending from -1 eV to 2 eV (in the minority spin 
channel these states are empty). At V both e g bands are de- 
generate, d 3z 2_ r 2 bands show a strong dispersion along T — Z, 
while d x 2_ y 2 bands are dispersive mainly along Z—A—R—Z. 
Both bands cross E-p at A and d 3z 2_ r 2 also halfway along 
R — Z. 



III. LAO/LNO SUPERLATTICES: ELECTRONIC 
PROPERTIES 

A. (LaA10 3 )i/(LaNi0 3 )i superlattices 

We now turn to the (LaA103)i/(LaNi03)i superlattices 
strained at the lateral lattice parameters of LAO, STO, and 
DSO. Their total density of states (DOS) are plotted in Fig. 3 
together with the one for bulk LNO. For compressive strain 
(short-dashed cyan line) the system is (semi-)metallic with a 
dip of DOS at the Fermi level. Moreover, due to the stack- 
ing of insulating LAO and conducting LNO layers an en- 
hanced resistivity is expected for (LaA103)i/(LaNi03)i at 
olao- In contrast, for tensile strain a gap in the total den- 
sity of states (Fig. 3) of 0.3 eV occurs above the Fermi 
level (solid line, gaey area). Interestnigly, the total DOS of 
(LaA103)i/(LaNi03)i at (Jdso almost coincides with the one 
for asTO around the Fermi level, implying a nonlinear depen- 
dence of the electronic properties on strain. Therefore, we 
concentrate in the following on the electronic properties of 
(LaA103)i/(LaNi03)i strained at olao an d osto- 
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TABLE I: Orbital occupation (total and majority/minority spin) of Ni e g states in (LAO) n /(LNO) n strained at asTO (tensile), cilao (compres- 
sive), and odso (tensile) compared to Ni in bulk LNO. The lateral strain with respect to the bulk LNO lattice parameter ao is e = (a — ao)/ao- 
Additionally, the orbital polarization of the occupied e g states P = (n x 2_ y 2 — n iz -i _ r i) f (n x n _ y i + n 3z 2_ r 2) and the magnetic moments 
are displayed. The last two columns give the number of holes Ap CT (difference in occupation with respect to the fully occupied orbitals) of 
the apical (a) and basal (b) oxygen within the MT sphere. 
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FIG. 3: (Color online) Total density of states for (LAO)i/(LNO)i for 
the case of tensile (STO) (solid line, gray area), (DSO) (long-dashed 
blue line) and compressive strain (LAO) (short-dashed cyan line). 
Additionally, the DOS of monoclinic bulk LNO is shown with a red 
line. The inset shows a zoom-in around the Fermi level displaying 
the opening of a band gap for (LAO)i/(LNO)i under tensile strain. 



In order to understand the origin of the distinct electronic 
behavior for tensile and compressive strain we analyze here 
the PDOS of the Ni 3d orbitals. Besides the filled t 2g orbitals 
there is a significant occupation of both e g orbitals at the Ni 
sites (see Fig. 4) split into bonding (between -7 and -5 eV) 
and antibonding contributions (from -1 to 2 eV). The most 
striking effect for tensile strain is a charge disproportionation 
on the Ni sublattice, Ni 3+ ^Ni 3+5 +Ni 3 ~' 531 where Nil has a 
stronger occupation of the e g orbitals. The effect is most pro- 
nounced for the majority spin states with 0.87e/0.86e (Nil) 
and 0.63e/0.71e (Ni2) in the d 3z 2_ r 2/d x 2_ y 2 orbitals, respec- 
tively. The reduced occupation in the majority spin channel 
at Ni2 is partially compensated by a stronger filling of mi- 
nority states at -4.3 eV. Altogether, the total difference in e g 
occupation at the Nil (Ni 3 ^) and Ni2 (Ni 3+(5 ) sites is — 0.2e, 
but there is a significant difference in the magnetic moments 
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FIG. 4: (Color online) Projected density of states (PDOS) of Ni 3d 
states (black dotted line), d 3z 2_ r 2 (green line) and d x 2_ y 2 orbitals 
(orange dashed line) in (LAO)i/(LNO)i for (a) compressive (LAO) 
and (b) tensile (STO) strain. The lower panels show the PDOS of 
basal (orange dashed line) and apical (green line) oxygens. Note the 
charge disproportionation at the Ni sites and the reduced density at 
the apical oxygen sites. 



1.45/ib (Nil) and 0.71/ib (Ni2) (for comparison, the mag- 
netic moment in LNO bulk is 1.12/zb). Without the octahe- 
dral tilts and rotations such a state is reported at U > 5 eV. 
A similar tendency for disproportionation was recently found 
in LNO films under tensile biaxial strain, 22 ' 23 but here as a 
result of confinement in the LAO/LNO superlattice the effect 
is strongly enhanced leading to the opening of a gap of ~0.5 
eV at the Fermi level. Concerning the orbital polarization, the 
two sites show opposite, largely compensating tendencies: a 
slightly higher occupation of d x 2_ y 2 at Ni2 (+2.2%) and the 
reverse effect at Nil (-1.7%). 

For compressive strain the charge disproportionation is no- 
ticeably reduced, reflected in a much smaller difference in the 
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magnetic moments (1.22 fi-g at Nil and 0.98 /ib at Ni2). We 
note that a non-disproportionated configuration lies only 12 
meV/20-atom unit cell higher in energy. The broadening of 
the conduction band leads to an overlap with the valence band 
which closes the band gap with a dip in the PDOS at Ep. As 
a consequence of the distortion (elongation) of the NiOg oc- 
tahedra discussed below, there is a stronger splitting between 
the d 3z 2_ r 2 and d x i_ y i states around Ep and an enhanced 
occupation of the d 3z 2_ r 2 orbital in the majority spin channel 
[0.86e (d 3z 2_ r 2) vs 0.76e (d x 2_ y 2) atNil and 0.77e (d 3z 2_ r 2) 
vs 0.72e (d x 2_ y 2) at Ni2]. Taking into account the occupa- 
tion of the minority spin chanel the total number of electrons 
within the MT sphere is 1.06 (d 3z 2_ r 2) vs. 0.98 (d x 2_ y i) at 
Nil and 1.02 (d3 Z 2_, r 2) vs 0.98 (d x 2_ y 2) at Ni2, respectively. 
This results in an orbital polarization of —3.8% and —1.6%. 
We note that much higher values were obtained by Han em- 
phet al. by integrating only over the antibonding states from 
1.5eVto^ F - 




(a) Majority spin channel Ni 1 
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FIG. 5: (Color online) Majority spin band structure of 
(LAO)i/(LNO)i under compressive strain (olao)- Note the reduced 
dispersion of the d 3z 2_,.2 band (green triangles) compared to bulk 
LNO caused by confinement and the opposite effect for the d x 2_ y 2 
band (orange circles) due to the reduction of the in-plane lattice pa- 
rameter. 

To gain more insight in the electronic structure of the 
(LaA103)i/(LaNi03)i superlattice, we have plotted in Fig. 5 
and Fig. 6 the band structure in the majority spin channel 
for compressive (olao) an d tensile (asTo) strain, respec- 
tively. As a consequence of confinement in the superlattice 
the d 3z 2_ r 2band has a reduced dispersion compared to bulk 
LaNiC>3 (Fig. 2) and the degeneracy of the two e g states at 
r is lifted with the e? 3z 2_ r 2 state being 0.3 (g,lao)/0.7 (osto) 
above Ep. The band structure around E-p is dominated by 
more localized e? 3z 2_ r 2 bands and strongly dispersive d x 2_ y 2 
bands (especially for compressive strain) with dispersion of 
the latter mainly along Z — A — R — Z. As already men- 
tioned in the discussion of the PDOS, there is a stronger split- 
ting between d 3z 2_ r 2 and d x 2_ y 2 states for olao compared to 
agio m contrast to bulk LNO (Fig. 2) a gap opens at A both 
for tensile and compresive strain. However, for compressive 
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FIG. 6: (Color online) Majority spin band structure 
(LAO)i/(LNO)i under tensile (STO) strain. d 3z 2_ r 2 and d x 2_ y 2 
character is emphasized by green triangles and orange circles, re- 
spectively. The direct band gap opens at A with a valence band max- 
imum (VBM) of predominantly Nil character and conduction band 
minimum (CBM) with mostly Ni2 contribution. 



strain (olao) the <i 3z 2_ r 2 band lies at the Fermi level along 
Z — A and halfway along R— Z causing a finite DOS at Ep. 
The upward shift of this band by ~ 0.5 eV opens a band gap 
for Osto- While the valence band maximum is mainly de- 
termined by d 3z 2_,.2 (and to a lesser extent d x 2_ y 2) states of 
Nil at A hybridized with p z orbitals of the apical oxygen, the 
conduction band minimum is dominated by d 3z 2_ r 2 states of 
Ni2, admixed with d 2 ,2_ J( 2and p z bands of the respective api- 
cal oxygen. This results in a stronger suppression of ligand 
holes of the apical oxygen of Nil. 

The charge disproportionation at the Ni sites for osto is 
visible in the spatial distribution of holes shown in Fig. 7. 
In contrast to the three-dimensional network of holes on the 
basal and apical oxygens in bulk LNO [Fig. 1(b)], a strik- 
ing feature in (LaA103)i/(LaNi03)i is the two- dimensional 
character of hole distribution within the LNO part of the su- 
perlattice due to a reduced number of holes in the p z orbital 
of the apical oxygen for both tensile and compressive strain 
(see also PDOS of oxygen in Fig. 4). At the same time, the 
in-plane oxygen orbital occupancies are very similar to those 
in bulk LNO. As a measure of the number of holes in Table I 
we have determined the difference in filling between p a and 
the fully occupied p n orbitals. We note that we consider here 
only the occupation within the MT sphere, which contains 60- 
70% of the total charge of oxygen. In bulk LNO there are 0.19 
less electrons in the p a than in the p^ orbital within the MT. 
This number is reduced to 0.09 — 0.15e for the apical oxygens 
in the superlattice. The suppresion is strongest for the api- 
cal oxygen of Nil (Ni 3 ~ d ) which simultaneously has a higher 
<i 3z 2_ r 2 occupation. Thus the suppression of holes due to the 
stronger ionicity of the Al 3+ -0-bond across the interface and 
the charge disporportionation at the Ni sites with more ionic 
and more covalent contributions at the two distinct Ni-sites 
are the crucial mechanisms in achieving an insulating state in 
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FIG. 7: (Color online) Electron density distribution of the unoccu- 
pied states integrated between Ef and Ef+2.5 eV in (LAO)i/(LNO)i 
for (a) compressive (LAO) and (b) tensile (STO) strain. Note the re- 
duction in hole density on the apical oxygen between Ni and Al, as 
well as the charge disproportionation at the Ni sites for tensile strain. 
La, O, and Al ions are shown in purple, red, and light gray, respec- 
tively. 
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FIG. 8: (Color online) (a) Projected density of states (PDOS) of 
3d states for Ni at interface (IF) and in the central layer (IF-1) in 
(LaA10 3 )3/(LaNi0 3 )3 for tensile (STO) strain. Total 3d, d 3z 2_ r 2 
and d x 2_ y 2 character are denoted by black dotted, green, and orange 
dashed line, respectfvelyl. (b) Electron density distribution of unoc- 
cupied states (yellow) integrated between Ef and Ef+2.5 eV. Note 
that only half of the simulation cell is shown. For color coding see 
Fig. 7. 



the interface, we observe again a reduced hole density at the 
apical oxygen at the interface. However, the 02p z hole den- 
sity is quickly restored to nearly bulk LNO behavior already 
in the next layer. This quick recovery of bulk LNO behavior 
at the apical oxygens away from the interface and the metallic 
state at the Ni sites with no evidence for charge disproportion- 
ation are instrumental for the insulator-to-metal transition in 
the (LNO)„/(LAO)„ SLs with increasing n, consistent with 
recent experiments. 4 ' 10 



B. (LaA103)3/(LaNi03>3 superlattices 

We turn now to the effect of LNO thickness on the elec- 
tronic properties of the superlattice. The projected density of 
states (PDOS) at the Ni sites in (LaA10 3 )3/(LaNi0 3 )3 at asTO 
is displayed in Fig. 8. Unlike the (LaA103)i/(LaNi03)i sys- 
tem, there is a finite DOS at Ep at the Ni sites both in the IF 
and IF-1 layer and no evidence for charge disproportionation. 
The e g states are nearly degenerate with a similar width as in 
the bulk and a small peak at E-p. Its character is predomi- 
nantly d 3z 2_ r 2 in the interface layer (IF) and Ni d x 2_ y 2 in IF- 
1 layer. As shown in the majority spin band structure in Fig. 9, 
these bands cross E-p at A [note that in (LaA103)i/(LaNi03)i 
a band gap opens at A]. The d 3z 2_ r 2 bands of the two distinct 
layers (IF and IF-1) are split, with the IF band being narrower 
and lying lower in energy. In contrast to recent soft x-ray re- 
flectivity measurements on (LaA103)4/(LaNi03)4 , we find 
only a small orbital polarization both in the interface (0.8%) 
and IF-1 (-0.3%) layers. 

Figure 8(b) shows the hole density distribution integrated 
between Ep and Ep+2.5 eV. Besides the similar hole density 
at the Ni sites independent of their position with respect to 
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FIG. 9: (Color online) Majority spin band structure of 
(LaA10 3 )3/(LaNi03)3 for tensile (STO) strain. Note the splitting of 
d 3z 2_ r 2 bands (green triangles) of Ni in the interface (IF) and (IF-1) 
layer, respectively. d x 2_ y 2 bands are emphasized by orange circles. 
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IV. STRAIN AND CONFINEMENT INDUCED LATTICE 
DEFORMATIONS 
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FIG. 10: (Color online) Strain dependence of the in-plane (red 
squares) and out-of-plane (blue circles) (a) bond angles and (b) Ni-O 
bond lengths in (LaA10a)i/(LaNi03)i; (c) Ni06 octahedral rotation 
(9) and tilt ((f)) angles, and (d) Glazer's rotation angles a and 7. The 
symbols at 0% strain represent the bulk LNO values. 

Not only the electronic but also the structural properties of 
the superlattice show a strongly asymmetric response to strain. 
The bond lengths and angles of the superlattices in compar- 
ison to bulk LNO and LAO are displayed in Table II and 
Fig. 10. For compressive strain the in-plane Ni-O bond length 
(1.92 A) in (LaA10 3 )i/(LaNi0 3 )i is significantly shorter than 
the out-of-plane distance (2.01 A), which is consistent with 
a volume-conserving tetragonal distortion and a crystal-field 
splitting between the d 3z 2_ r 2 and d x 2_ y 2 orbitals observed in 
Fig. 4(a). In contrast, the charge disproportionation for ten- 
sile strain goes hand in hand with bond disproportionation 
driven by a breathing mode with two nearly regular octahe- 
dra with in-/out-of -plane cation-anion distances of 2.02/1.99 
A and 1.93/1.93 A at the Nil and Ni2 sites at asTO, respec- 
tively. 

As mentioned previously, in the bulk both LAO and LNO 
attain a rhombohedral structure with rotations of the octahe- 
dra around all three pseudocubic axes by a = ft = 7 = 3.04° 
(LAO) and 5.4° (LNO) . This corresponds to Glazer's no- 
tation (a~, a - , a~). Recent DFT calculations have found 
that under compressive strain the in-plane rotations of LAO 
are suppressed and the rotation pattern changes to (a , a , c~). 
Vice versa for tensile strain the rotations are around the [110] 
direction (a~ ,a~ ,c°). In contrast, an (a~ ,a~ ,c~) pattern 



was reported in a combined XRD and DFT study on strained 
LNO films with a strong reduction of 7 = 0.3 ± 0.7° and 
enhanced a = 7.1 ± 0.2° for tensile (agTo) and an oppo- 
site trend for compressive ((Jlao) strain [a = 1.7 ± 0.2° and 
7 = 7.9 ± 0.9°). While the general pattern of octahedral ro- 
tations is preserved in the (LaA103)i/(LaNi03)i superlattice, 
the presence of the interface and the different strain state in 
the LAO and LNO part imposed by the substrate lead to sig- 
nificant changes of the absolute values of rotation angles com- 
pared to the individual phases described above. The rotational 
angle of the NiOg octahedra around [100] increases with strain 
from a ~ 4.7° (a LAO ) to ~ 6.1° (a S To) and ~ 7.2° (o D so), 
whereas 7 remains constant ~ 6.2°. On the other hand, for 
the AIO6 octahedra the rotational angle around [001] changes 
from 7 = 4.4° (compressive) to 7 ~ 2.5° (tensile strain). 

The cation-oxygen-cation bond angles have a direct influ- 
ence on the band width. In bulk LNO the Ni-O-Ni angle is 
164.8°. In LAO, the deviation of the Al-O-Al angle from 
180° is smaller (171.4°). For compressive strain: the out- 
of-plane Ni-O-Al angle in (LaA103)i/(LaNi03)i is increased 
to ~ 166 — 167°, while the in-plane Ni-O-Ni angle decreases 
to 163° (the in-plane Al-O-Al angle, 165.6°, is also reduced 
w.r.t. the bulk value). In contrast, for tensile strain both the 
in- (Ni-O-Ni) and out-of-plane (Ni-O-Al) bond angles are sig- 
nificantly reduced to ~ 162° for osto and ~ 161 — 162° 
for Odso- The Al-O-Al is 167.3°, again smaller than the 
corresponding angle in bulk LAO (171.4°). The decrease of 
bond angles reduces the band width of d 3z 2_ r 2 and d x 2_ y 2 
bands between E-p and Ep+4 eV, as can be seen in the PDOS 
in Fig. 4b). The Ni-O-Ni angle has been previously identi- 
fied as the control parameter of the MIT in rare earth nicke- 
lates, where a smaller deviation from 180° favors Ni 3d-02p 
hybridization and a metallic state, while a strong deviation 
leads to a MIT . Indeed, the values we obtain for tensile 
strain (160 - 162°) are close to the Ni-O-Ni angles in PrNi0 3 
(158.7°) and NdNiOg (157.1°) 25 that undergo a MIT. 

In the (LaA103)3/(LaNi03)3 superlattice under tensile 
strain (agTo) the octahedra are nearly regular with an in-/out- 
of-plane Ni-O bond length of 1.97/1.95 A (IF) and 1.97/197 A 
(IF-1), respectively. The pattern of octahedral rotations is pre- 
served with almost identical absolute values of rotation angles 
compared to the (LaA103)i/(LaNi03)i at asTO The out-of- 
plane Ni-O-Al (164°) and Ni-O-Ni (162°) and in-plane Ni-O- 
Ni angles - 162° (IF)/163° (IF-1) are slightly higher than the 
ones for the (LaA103)i/(LaNi03)i superlattice in line with 
the metallicity of the system. 



V. SUMMARY 

While in the bulk LNO is metallic at all temperatures with 
degenerate e g orbitals, we find that quantum confinement to- 
gether with tensile strain allows to stabilize an insulating state 
in a (LaA103)i/(LaNi03)i superlattice. The charge dispro- 
portionation at the Ni sublattice with one strongly and one 
weakly magnetic Ni ion is fostered by the lattice response to 
strain, expressed in decreasing bond angles and enhanced tilt 
and rotation angles that reduce the band width. The charge 
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TABLE II: Ni-oxygen out-of-plane (d,Ni-o a ) and in-plane (d Ni _ b) distances, out-of-plane bond angle 9 C (Al-O-Ni), and in-plane angle 6 a b 
(Ni-O-Ni). Additionally, the rotation 9 and tilt cj> angles, as well as Glazer's rotation angles a and 7, are displayed. 







dNi-O a 






Sab 


6 





a 


7 


l/l@a L AO 


Nil 


2.02 


1.93 


166.9 


162.5 


6.2 


6.5 


4.6 


6.2 


e =-1.1% 


Ni2 


2.00 


1.91 


166.4 


162.5 


6.2 


6.8 


4.8 


6.2 


l/l@asTO 


Nil 


1.99 


2.02 


163.7 


162.4 


6.1 


8.1 


5.8 


6.1 


e =1.7% 


Ni2 


1.93 


1.93 


162.0 


162.1 


6.1 


9.0 


6.4 


6.1 


l/l@aDSO 


Nil 


2.00 


2.05 


160.7 


161.7 


6.2 


9.7 


6.8 


6.2 


e =2.7% 


Ni2 


1.94 


1.95 


158.3 


161.0 


6.2 


10.8 


7.7 


6.2 


3/3@asTO 


Ni IF 


1.95 


1.97 


163.9 


162.3 


6.2 


8.1 


5.7 


6.2 


e =1.7% 


NiIF-1 


1.97 


1.97 


162.4 


162.6 


6.4 


8.8 


6.2 


6.4 


LaNiOs 


Ni 


1.94 


1.94 


164.8 


164.8 


5.4 


7.6 


5.4 


5.4 


LaAlOs 


Al 


1.90 


1.90 


171.4 


171.4 


3.0 


4.3 


3.0 


3.0 



ordered state in the (LaA103)i/(LaNi03)i SL under tensile 
strain is characterized by two nearly regular octahedra with 
distinct bond lengths and hence a reduced crystal field split- 
ting at the Ni-sites. On the other hand, under compres- 
sive strain bands broaden and close the band gap and CD is 
strongly inhibited. Expansion of the out-of-plane Ni-O bond 
length results in a crystal-field splitting between the d x 2_ y 2 
and d 3z 2_ r 2 orbitals and an energy lowering of the latter [see 
Fig. 4(b)]. This is similar to previous results for LNO films 
under compressive strain. 

LAO/LNO superlattices and strained LNO films show simi- 
lar tilt patterns but different magnitude of octahedral rotations. 
This gives further evidence that the lattice response to strain 
in superlattices is very different from the one in thin films and 
thus provides a further possibility to tune the band width and 
the electronic behavior in transition metal oxides, as recently 
reported for LNO/SMO SLs. 28 

The large band width of the e g states precludes signifi- 
cant changes in the orbital occupation via strain. The lack 
of lifting of orbital degeneracy is similar to other nickelates, 
e.g., NdNiOa where the MIT is not associated with orbital 
ordering. 29 ' 30 Recently, Mazin emphet al. proposed that the 



CO state in rare-earth nickelates gets gradually suppressed 
under pressure due to overlap of bands. In the confined 
(LaA103)i/(LaNi03)i we observe an analogous "melting" of 
the CD state as a function of strain: The significant difference 
in magnetic moments at the two distinct sites Nil and Ni2 
sites at osto is gradually reduced for compressive strain. 

The insulating state in the (LaA103)i/(LaNi03)i super- 
lattice is achieved through a combination of charge dispro- 
portionation and a reduction of the hole density on the api- 
cal oxygens at the interface due to the ionic bonding to Al. 
The metallic behavior at the Ni sites and the hole density 
of the apical oxygen in inner layers is quickly restored in 
(LaA103)3/(LaNi03)3 superlattices, in line with recent obser- 
vations of a MIT with increasing LNO thickness. • We note 
that due to the presence of a surface, the mechanism of MIT 
in single layer films • may be different and requires further 
investigation. 
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